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ABSTRACT 

Magnetospheres of pulsars are thought to be filled with plasma, and variations in plasma supply 
can affect both pulsar emission properties and spin-down rates. A number of recently discovered 
"intermittent" pulsars switch between two distinct states: an "on", radio-loud state, and an "off", 
radio-quiet state. Spin-down rates in the two states differ by a large factor, ~ 1.5 — 2.5, which is not 
easily understood in the context of current models. In this Letter we present self-consistent numerical 
solutions of "on" and "off" states of intermittent pulsar magnetospheres. We model the "on" state 
as a nearly ideal force-free magnetosphere with abundant magnetospheric plasma supply. The lack 
of radio emission in the "off" state is associated with plasma supply disruption that results in lower 
plasma density on the open field lines. We model the "off" state using nearly vacuum conditions on 
the open field lines and nearly ideal force- free conditions on the closed field lines, where plasma can 
remain trapped even in the absence of pair production. The toroidal advection of plasma in the closed 
zone in the "off" state causes spin-downs that are a factor of ~ 2 higher than vacuum values, and we 
naturally obtain a range of spin-down ratios between the "on" and "off" states, ^ 1.2 — 2.9, which 
corresponds to a likely range of pulsar inclination angles of 30—90°. We consider the implications 
of our model to a number of poorly understood but possibly related pulsar phenomena, including 
nulling, timing noise, and rotating radio transients. 

Subject headings: magnetohydrodynamics — pulsars: general — stars: magnetic field 



1. INTRODUCTION 

Pulsars spin down due to torques exerted by currents 
flowing on the surface of the neutron star. In the absence 
of magnetospheric plasma, t he pulsar sp i ns down due 
tomagneto-dipole radiation (Michel 199l| iBeskin et al.l 
119931 ). Magnetospheric plasma allows for c urrents that 
can p roduce additional spin-down torques (jSpitkovskvl 
120061 hereafter S06), so variations in plasma supply can 
potentially modulate pulsar spin-down. 

Several classes of pulsars have been identified that may 
exhibit such modulation. Nulling pulsars have radio 
emission that appears to shut off for a few to several 
tens of rotation periods. Since radio emission is pre- 
sumably tied to the magnetospheric plasma, the nulling 
suggests that some process is affecting the plasma supply 
and/or currents above the polar ca ps (|Wang et al.ll2007t 
iZhang et al]l2QQ7t lTimokhinll2010l ). Intermittent pulsars 
switch between an "on" , radio-loud, state in which they 
behave like normal radio pulsars, and an "off", radio- 
quiet, state in which they produce no detectable radio 
emission for long periods of time. This process may be 
an extreme manifestation of nulling. The first two inter- 
mittent pulsars with publishe d data have quite di fferent 
duty cycles: PSR B1931+24 (jKramer et al.ll2006l here- 
after K06), with a period P « 0.8 s, cycled through 
the "on" -"off" sequence of states a pprox i matel y once a 
month, whereas PSR J1832+0029 (|Lvnd I2009D . with a 
period P w 0.5 s, kept quiet for nearly two years be- 
tween "on" cycles of unknown length. The spin-down 
rate for each of these pulsars is larger in the "on" state 
than in the "off" state by a factor /on^off — 1-5. A third 
intermittent pulsar, PSR J1841-0500, was not detected 
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for over 1.5 years between "on" cycles, on e of which ap- 
pears to have lasted for at least a year (|Camilo et al.l 
120111 ). This pulsar seems to have a spin-down ratio be- 
tween "on" and "off" states of / on ^off — 2.5. Such sub- 
stantial differences in spin-down rates suggest that the 
pulsar magnetosphere undergoes a dramatic reconfigura- 
tion as it transitions between the "on" and "off" states, 
yet such a transition was reported for PSR B1931+24 to 
take place in just over 10 pulsar periods. 

Intermittent pulsars offer a unique testbed of pulsar 
theory. K06 first proposed that in the "on" state plasma 
fills the pulsar magnetosphere and supports plasma pro- 
cesses that produce radio emission. The pulsar transi- 
tions to the "off" state when open field lines become 
depleted of charged radiating particles. K06 approx- 
imated the spin-down rate in the "off" state by the 
spin-down of a vacuum dipole and estimated the ex- 
tra plasma currents needed to account for the observed 
spin-down of the "on" state. It is not clear, however, 
whether a working pulsar should naturally yield the re- 
quired on-off spin-down ratio of this picture. The sim- 
plest model for the "on" state is the force-free magne- 
tosphere, which has abundant charges everywhere. The 
force-free spin-down rate is larger than the vacuum spin- 
down rate by a factor /ff^ V ac = (1 + sin 2 a)/ (2/3 sin 2 a) 
(S06) that is greater than or equal to 3 for all inclination 
angles a. This is clearly inco mpatible with the observed 
values, .forwoff ~ 1.5 - 2.5 (|Beskin fc Nokhrinal [20071 : 
iGurevich fc Istominl 12007). This suggests that, perhaps, 
we do not understand the spin-down power of the "off" 
state. 

The shutoff of pair formation as the intermittent 
pulsar switches "off" allows plasma to escape along 
the open field lines, but the plasma in the closed 
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zone is confined by the geometry of the field lines. 
This is an important physical effect that w as not in- 
cluded in previous work (KQ6: iLi et alJ [20TTl hereafter 
LSTll lKalapotharakos etaH l2011h . The currents and 
charges associated with plasma trapped in the closed 
zone can increase the spin-down in the "off" state even 
if the open field lines are empty. Thus, in this paper, we 
model the "off" state using a simple two-zone prescrip- 
tion in which the closed zone is highly conducting and the 
open field lines are vacuum-like. We run resistive force- 
free simulations to directly solve for the magnetospheric 
geometry in the "on" and "off" states and test whether 
this model can produce the observed intermittent pulsar 
spin-down ratios. In Section [2] we describe the numeri- 
cal code and setup. Section [3] illustrates our intermittent 
pulsar solutions and shows the spin-down results. Sec- 
tion 3] provides a brief summary of our results and their 
observational implications. 

2. SETUP 

We employ a three-dimensional numerical code (see 
S06) that imple ments the finite differenc e time-domain 
scheme (FDTD, iTaflove & Hag ness 2005) to evolve elec- 
tromagnetic fields from Maxwell's equations, 

d t E = cV x B - 4ttJ, 

d t B = -cV x E, (1) 

where the current j is given by 

j = pv + aEft uid . (2) 

The fluid velocity v = c(E x B)/(B 2 + Eq) is the general- 
ized drift velocity, £Wid = ^(E + vxB) is the fluid frame 
electric field, 7 = (1 — v 2 /c 2 ) -1 / 2 , Eq is the magnitude 

of i^fluid (see LST11), p is the charge density, and a is 
the plasma conductivity in the fluid frame. The central 
region of our grid is occupied by a conducting spherical 
star of radius i?*, rotating at angular velocity f2, with 
embedded dipole field of magnetic moment ft inclined 
relative to the rotation axis by angle a. We resolve the 
light cylinder R^c = c/fl with 80 cells and set R* = 30 
cells. See LST11 for a detailed description of our code 
and resistive current formulation. We have verified that 
our solutions are converged with spatial resolution, as 
well as run sufficiently long so as to reach a steady state 
in the frame corotating with the pulsar. 

In LST11 we showed that our resistive force- free formu- 
lation can capture both the vacuum and ideal force-free 
limits by varying the conductivity parameter. We model 
the "on" state as a magnetosphere with high conductivity 
(<j/Q) 2 = 40, our fiducial value representing force- free- 
like conditions for i?*/i?LC = 3/8. This is preferable to 
using an ideal force- free formulation as in S06, because 
our high conductivity solutions are numerically cleaner, 
especially in the current sheets (LST11). We model the 
"off" state as having conducting and vacuum-like regions 
separated sharply at the boundary of the closed field line 
region, which we approximate as the closed field lines 
of the force- free "on" state. Nonrotating dipole mag- 
netic field lines ar e traced by the curves s = sq sin 2 
(jMichel fc Lil Il999h . where is the angle from the mag- 
netic axis, s is spherical radius, and sq> the maximum 



perpendicular distance of a field line from the magnetic 
axis, specifies which field line is under consideration. The 
closed field lines in the force-free solutions are typically 
stretched in the direction perpendicular to the magnetic 
axis as compared to vacuum dipole closed field lines. We 
use the surface s = R\,c sin 2 #(i?LC sin 2 / R*) k to de- 
marcate the closed field line region. The exponent k is 
set at each inclination angle to best match the shape of 
the force- free closed zone (see Section [3] for illustrations 
of these demarcation surfaces), and the half-angle size of 
the conducting polar cap is given by pc = y/R*/ Rlc • 
Interior to the demarcation surface, we set (cr/Q) 2 = 40, 
as in the "on" state. Exterior to the demarcation sur- 
face, we set (cr/Q) 2 = 0.04, a fiducial value representative 
of vacuum-like conditions for R^/R^c =3/8. We pick 
this value over <j/Q = for better numerical accuracy 
when computing the field geometry. This choice of con- 
ductivity on open field lines and the exact shape of the 
demarcation surface have minimal effect on spin-down 
(see Section [3]). 

3. RESULTS 

Figure [T] shows magnetic field lines in the jl — Cl plane 
for inclined dipoles at a = 30° (top), 60° (middle), 90° 
(bottom). Color is representative of the out-of-plane 
magnetic field (as in LST11). The left column shows 
the "off" state with abundant plasma in the closed zone 
but a shortage of plasma along open field lines. The red 
curves indicate the boundary between the conducting 
and vacuum-like regions in the illustrated cross-section 
of the magnetosphere. The right column shows the "on" 
state with abundant plasma everywhere. The magneto- 
spheres in the "on" state are force- free-like, with conduc- 
tion currents flowing along open field lines. The current 
returns through the current sheets and along the bound- 
ary of the closed field line region in the current layers. 
Gross features of the magnetosphere in the "off" state are 
vacuum-like, with large closed field line region and dis- 
placement currents contributing to spin-down. For refer- 
ence we show a representative vacuum solution with in- 
clination angle a = 60° in Figure [21 There are a number 
of important differences that distinguish the intermittent 
"off" state from the vacuum solutions. The current from 
toroidal advection of charged plasma in the closed zone 
leads to greater magnetic flux passing through the light 
cylinder and a larger fraction of open field lines. Fur- 
ther, poloidal conduction currents are present even out- 
side the conducting closed zone. They are due primarily 
to the fluid advection term pv in the current and lead 
to greater magnetic field sweepback and stronger current 
sheets than in pure vacuum solutions. 

Figure [3] shows the spin-down luminosity, normalized 
by Lq (where Lq is 3/2 times the power of the orthogonal 
vacuum rotator with finite R*/Rlc — 3/8), as a function 
of inclination angle for both the "on" and "off" states of 
the magnetosphere. Spin-down luminosity is calculated 
as the surface integral of the Poynting flux over a sphere 
of radius i?LC • We also show for reference the spin-down 
for the vacuum solution in the total absence of plasma. 
The "on" state spin-down is essentially the force-free 
spin-down. There is an uncertainty of roughly 10% of 
Lq in the force-free spin-down at low inclination angles 
due to stellar boundary effects and unphysical dissipa- 
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Fig. 1. — Magnetic field lines in the ft — f2 plane for inclined dipoles at a = 30° (top), a = 60° (middle), a = 90° (bottom). Color is 
representative of out-of-plane magnetic field into (red) and out of (blue) the page. The left column shows the "off" state with abundant 

plasma in the closed zone and vacuum-like conditions along open field lines. The red curves indicate the intersection with the ft — Q plane of 
the demarcation surface between these distinct conducting and vacuum-like regions. Gross features of the magnetosphere in the "off" state 
are vacuum-like, but the "off" states have more magnetic flux passing through the light cylinder and stronger magnetic field sweepback 
than vacuum solutions. The right column represents the "on" state with abundant plasma everywhere. 



4 




Fig. 2. — Magnetic field lines in the jl — Q, plane for a vacuum 
inclined dipole at a = 60°. Color is representative of out-of-plane 
magnetic field into (red) and out of (blue) the page. There is large 
closed field line region, and displacement currents circulate and 
contribute to spin-down. 



tion above the polar caps (see LST11). The spin-down 
in the "off" state lies between the force-free and vac- 
uum spin-down values at all inclination angles. The spin- 
down of the aligned rotator in the "off" state is small, as 
the displacement currents are zero and large-scale con- 
duction currents are weak. The open field lines carry 
minimal Poynting flux. Inclined dipole open field lines 
carry Poynting flux, however, and the larger magnetic 
flux passing through the light cylinder in the "off" state 
results in higher spin-down than in the vacuum solution. 
This increase in Poynting flux over vacuum spin-down 
in the "off" state is largest at high inclination angles, 
where the open field lines in the vacuum solutions carry 
the most Poynting flux. The advective poloidal conduc- 
tion currents also lead to higher spin-down since they 
cause larger magnetic field sweepback than in vacuum so- 
lutions. We evaluated the relative contributions of these 
two effects to higher spin-down over vacuum solutions by 
killing the advection current term pv outside the closed 
zone. The larger magnetic flux passing through the light 
cylinder is responsible for upwards of two-thirds of the 
increase in Poynting flux over vacuum spin-down in the 
"off" state. 

The spin-down in the simulated "off" state is relatively 
insensitive to the size of the conducting closed zone. We 
varied the volume of the conducting closed zone by a 
factor of a few while keeping the size of the polar cap 
fixed, and we find that spin-down changes by less than 
20% of Lo at all inclination angles. It is not the total 
size of the conducting closed zone, but rather its extent 
on the stellar surface that is more important in deter- 
mining the gross magnetospheric properties. Blocking 
radial currents near the stellar surface limits the large- 
scale poloidal current circuit. Spin-down also depends 
very weakly on the fiducial values of a/Cl that we pick 
to represent highly conducting and vacuum conditions. 
We checked this by setting the conductivity to a/Q = 
instead of (cr/Q) 2 = 0.04 along open field lines, and we 
find that the resulting spin-down is offset by less than 
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Fig. 3. — Spin-down luminosity dependence on inclination angle 
for the force-free-like "on" state, the "off" state, and the vacuum 
solution. The "off" state spin-down naturally lies between the 
force- free and vacuum spin-down values for all inclination angles, 
spin-down is normalized by 3/2 times the spin-down power of the 
orthogonal vacuum rotator. 



20% of Lo for all inclination angles. The aligned "off" 
state spin-down drops to zero in this case, as we expect, 
since the large-scale conduction currents have been elim- 
inated. We also tried setting the closed zone to have 
conductivity (cf/Q) 2 = 20 instead of (cf/Q) 2 = 40, and 
spin-down results change by less than 10% of Lo at all 
inclination angles. 

The most relevant observational parameter in the con- 
text of intermittent pulsars is the ratio of spin-down 
power in the "on" and "off" states, / on -H>ff- Figure |4] 
shows this ratio for our models. We emphasize that 
these models provide a physically well-motivated set of 
solutions to describe both the "on" and "off" states of 
intermittent pulsars. Assuming a uniform distribution 
of pulsar inclination angles, we expect two- thirds of in- 
termittent pulsars to have f n->off between 1.2 — 2.9. 
Known intermittent pulsars fa ll within this range (K06, 
lLvndl2009l ; iCamilo et al.ll2011f ). Uncertainties in the ex- 
act charge configuration and charge transport properties 
in the "off" state magnetosphere imply error bars as- 
sociated with our calculation of the ratio / n-H>ff- The 
uncertainty is of order 10% at inclination angle a = 90° 
and rises with decreasing inclination angle to a factor of 
order 2 at a = 30°. Below a = 30° the ratio / on ^ofT rises 
above 3 due to the small spin-down of the aligned rotator 
in the "off" state, but there are large uncertainties in the 
exact value of / n-K>ff i n this regime. 

A major improvement of our intermittent pulsar mod- 
els over existing work (LST11) lies in the treatment of 
accelerating potential drops, which can be used as a fidu- 
cial measure of energy gain as particles fly away from 
the stellar surface. We define the potential drops as in 
LST11, i.e., as the line integral along magnetic field lines 
of the corotating electric field E' = E + (f2 x r) x B / c. In 
ideal force-free solutions the potential drops vanish, as 



Intermittent Pulsars 



5 



20 40 60 80 

Inclination angle, degrees 

Fig. 4. — The ratio of spin-down power between the "on" and 
"off" states, /o n ^off ; as a function of inclination angle. We obtain 
/on^off ~ 1-2 — 2.9 for a > 30°. Known intermittent pulsars have 
/on-^off in this range. 



E ■ B = 0, but they increase monotonically with increas- 
ing bulk resistivity. Previously, we implemented a con- 
stant conductivity a/Vt throughout the magnetosphere. 
As the polar cap played no fundamental role in these 
models, open field line potential drops were generally 
limited by the full pole-to-equator potential drop, yield- 
ing unphysically large potential drops of order 10 16 V. 
Our new models for the "off" state introduce conducting 
plasma in the closed field line zone, effectively shielding 
the potential drop there. The potential drops are then 
limited by the polar cap potential drop V pc = |/?|/-Rlc- 
For typical pulsars with periods P ~ Is, we obtain char- 
acteristic potential drop of ~ 10 12 V, more in line with 
expectations. 

4. DISCUSSION 

We have developed an improved numerical method 
for solving pulsar magnetospheres with resistivity, and 
we apply it to describe intermittent pulsars in a self- 
consistent manner. In the "on" state, plasma is abundant 
everywhere, and the instabilities in this plasma gener- 
ate coherent radio emission. In the "off" state plasma 
has leaked off the open field lines, suppressing the open 
field line currents. The radio emission is hence shut off. 
Plasma remains trapped in the closed zone, however, and 
the current from toroidal advection of these charges leads 
to spin-down values that are a factor of ~ 2 larger than 
vacuum values. This allows us to naturally produce spin- 
down ratios / on ^off ^ 1-2 — 2.9 for inclination angle in 
the range a = 30° — 90°, consistent with observations, 
and we obtain realistic values for accelerating potential 
drops, Vdrop ~ 10 12 G. In our model the spin-down ratio 
/on-^off takes on its minimum value at a = 90° and in- 
creases monotonically with decreasing inclination angle. 
Hence, given an observed spin-down ratio / on -^off f° r an 
intermittent pulsar, we can predict, within errors, the 
pulsar's inclination angle. Alternatively, given the pul- 



sar inclination angle from, e.g., polarization vector sweep 
data, we can predict the spin-down ratio f on ->off if the 
pulsar displays clear "on" and "off" states. A verification 
of our predictions would lend strong support to the idea 
that abrupt changes in plasma supply on open field lines 
can play an important role in determining the emission 
properties and spin-down rates of real pulsars. 

One potential limitation of our models relates to 
the spatial distribution of conductivity in the mag- 
netosphere. We specify only a conducting torus lo- 
cated in the plane perpendicular to the magnetic axis. 
In fact, rotation of a conducting neutron star leads 
to unipolar induction, and domes of negative charge 
form above the magnetic poles even in the absence of 
a pa ir cascade if the work function of the surface is 
lo^jKraaise^Polstorff & Mich el 11984 119851: IPetrietalJ 
l2002aI FSpitkov skv fc AromJ 12002b ISpitkovskvl 12004( 1. 
We explored this effect by prescribing additional con- 
ducting domes above the magnetic poles, setting 
(<j/Q) 2 = 40 interior to the surfaces specified by s = 
i?LC cos 2 6. Spin-down values increase at all inclination 
angles by an offset of less than 20% of Lq, but the spin- 
down ratio /on^off remains within the errors quoted in 
Section 03 

It is important to note that at present the precise 
mechanism by which the plasma supply is turned on and 
off is still unclear. In at least one intermittent pulsar 
case the switching between different spin-down rates is 
quasi-periodic (K06), implying that the processes supply- 
ing and limiting plasma alternately recur. Inter mittency 
may actually be the same basic process as nulling, the 
crucial difference being that the "off" state lasts months 
to years instead of from a few rotation periods up to 
days. In this light, timing noise can in some instances 
be caused by nulling events that are not resolved by the 
observations. The pulsar timing data in the Jodrell Bank 
data archive are typically smoothed and do not resolve 
features that occ ur on time scales s horter than a few 
tens of days (see iHobbs et al.l 120101 : iLvne et all 12010( 1 . 
Suppose recurring nulling events last for an accumulated 
time that is of order a few percent of the time that the 
pulsar behaves normally. Since in our model the spin- 
down in the "off" state is typically of order 1/2 that 
in the "on" state, the resulting spin-down luminosity in 
a state wi th unresolved mill ings could easily be modified 
by ~ 1%. ILvne et "all (|2010l ) established a connection be- 
tween timing noise and mode changing, when the pulse 
profile of the pulsar changes to a different shape, but 
they do not see mode changing in all pulsars that exhibit 
timing noise. It is possible that some of these pulsars 
undergo nulling events, which leads to variation in spin- 
down rate and the observed timing noise. The mode 
changing is presumably due to changes in the magne- 
tospheric configuration. One possibility is that it is an 
intermediate state between our "on" and "off" states in 
which some but not all of the plasma along open field 
lines has leaked away. These intermediate states could 
lead to a number of observed pulsar features including 
jitter, subpulse drift, precursors, and interpulses, though 
such ideas are at present still speculative. Another pos- 
sible cause for changing pulse profiles is the presence of 
multiple components to the emission, e.g., the core and 
conal components, and one or more components being 
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modified or shut off by changes in plasma supply. 

If the "on" and "off" states of the pulsar can last any- 
where from a few periods to years, it is further con- 
ceivable that Rotating Radio Transients (RRATs) are 
another manifestation of the process captured by our 
model. RRATs are rotating neutron stars that occa- 
sionally emit pulses, usually is olated, but in a few in- 
stances in a string of several ([McLaughlin et all 120061 : 
iPallivaguru et af]l2011f ). It has not been ruled out that 
RRATs are an extreme form of nulling, and we envi- 
sion RRATs as pulsars that spend the majority of their 
time in our "off" state, but turn "on" from time to time 



and emit pulses. If future observations are able to es- 
tablish a more definitive link between the seemingly dis- 
parate processes of intermittency, nulling, timing noise, 
and RRATs, it would be an important step in our un- 
derstanding of pulsar physics. 
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ing Center and Visualization Laboratory. 



REFERENCES 



Beskin, V. S., Gurevich A. V., Istomin, Ya. N. 1993, Physics of 
the Pulsar Magnetosphere (Cambridge: Cambridge University 
Press) 

Beskin, V. S. & Nokhrina, E. E. 2007, Ap&SS, 308, 569 
Camilo, F., Ransom, S. M., Chatteriee, S., Johnston, S., Sz 

Demorest, P. 2011, arXiv: 1111.5870 
Gurevich, A. V., Sz Istomin, Y. N. 2007, MNRAS, 377, 1663 
Hobbs, G., Lyne, A. G., & Kramer, M. 2010, MNRAS, 402, 1027 
Kalapotharakos, C, Kazanas, D., Harding, A., & Contopoulos, I. 

2011. laTXivTl l08.2138 
Kramer, M., Lyne, A.G., O'Brien, J.T., Jordan, C.A., Lorimer, 

D.R. 2006, Sci, 312, 549 
Krause-Polstorff, J., & Michel, F. C. 1984, MNRAS, 213, 43 

Krause-Polstorff, J., & Michel, F. C. 1985, A&A, 144,72 

Li, J., Spitkovsky, A., & Tchekhovskoy, A. 2011, arXiv: 1107.0979 
Lyne, A. G. 2009, ASSL, 357, 67 

Lyne, A., Hobbs, G., Kramer, M., Stairs, I., & Stappers, B. 2010, 

Science, 329, 408 
McLaughlin, M. A., Lyne, A. G., Lorimer, D. R., et al. Nature, 

439, 817 

Michel, F. 1991, Theory of Neutron Star Magnetosphere 
(Chicago: University of Chicago Press) 



Michel, F. C. & Li, H. 1999, Physics Reports, 318, 227 
Palliyaguru, N. T., McLaughlin, M. A., Keane, E. F., et al. 2011, 
MNRAS, 1871 

Petri, J., Heyvaerts, J., & Bonazzola, S. 2002a, A&A, 384, 414 
Petri, J., Heyvaerts, J., & Bonazzola, S. 2002b, A&A, 387, 520 
Spitkovsky, A. 2004, IAUS, 218, 357 
Spitkovsky, A. 2006, ApJ, 648, L51 

Spitkovsky, A., Arons, J. 2002, in 'Neutron Stars in Supernova 

Remnants', eds P. O. Slane and B. M. Gaensler, (San 

Francisco: ASP), 271, 81 
Taflove, A. N. & Hagness, S. C. 2005, Computational 

Electrodynamics: The Finite-Difference Time-Domain Method, 

Third Edition (Norwood, MA: Artech House) 
Timokhin, A. N. 2010, MNRAS, 408, 41 

Wang, N., Manchester, R. N., & Johnston, S. 2007, MNRAS, 377, 
1383 

Zhang, B., Gil, J., & Dyks, J. 2007, MNRAS, 374, 1103 



